Prepared for submission to JCAP 



Enhanced Preheating after Multi-Field 
Inflation: On the Importance of being 
Special 

(N 

o 

Thorsten Battefeld 1 ) Alexander Eggemeier 1 ) John T. Giblin, Jr. 2 )' 3 ) 

in 

1) Institute for Astrophysics, University of Goettingen, Friedrich Hund Platz 1, D-37077 
y—i Goettingen, Germany 

2) Department of Physics, Kenyon College, Gambier, OH 43022, U.S.A. 

3) Department of Physics, Case Western Reserve University, Cleveland, OH 44106, U.S.A. 

E-mail: tbattefe@astro.physik.uni-goettingen.de, 
a.eggemeier@stud.uni-goettingen.de, giblinj@kenyon.edu 

Qh 

Abstract. We discuss preheating after multi-field inflation in the presence of several preheat 
+3 matter fields that become light in the vicinity of (but not at) the inflatons' VEV, at dis- 

tinct extra-species-points (ESP); this setup is motivated by inflationary models that include 
particle production during inflation, e.g. trapped inflation, grazing ESP encounters or mod- 
ulated trapping, among others. While de-phasing of inflatons tends to suppress parametric 
resonance, we find two new effects leading to efficient preheating: particle production during 
the first in-fall (efficient if many preheat matter fields are present) and a subsequent (nar- 
row) resonance phase (efficient if an ESP happens to be at one of several distinct distances 
from the inflatons' VEV). Particles produced during the first in-fall are comprised of many 
ON species with low occupation number, while the latter are made up of a few species with high 

occupation number. 

We provide analytic descriptions of both phases in the absence of back-reaction, which 
we test numerically. We further perform lattice simulations to investigate the effects of back- 
reaction. We find resonances to be robust and the most likely cause of inflaton decay in 
r> multi-field trapped inflation if ESP distributions are dense. 
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1 Introduction 

After inflation the Universe is empty and cold, in dire need to be reheated so that primordial 
nucleosynthesis can take place. The decay of the inflaton(s) can be perturbative (reheating) 
or explosive (preheating), depending on the inflationary model (see [1-3] for reviews). 

Preheating after inflation is often treated phenomenologically, since the degrees of free- 
dom to which the inflaton(s) couple are unknown. Exceptions are certain low scale inflation- 
ary models, such as reheating after inflation within the MSSM [4] or Higgs inflation [5-7]. 
It was soon discovered that non-perturbative effects, such as parametric resonance [8-10] or 
tachyonic instabilities [11, 12], can lead to a fast transfer of energy from the inflaton sector 
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to a preheat matter field. These effects, if present, can dominate over the decays channels 
of the old theory of reheating [13, 14] 1 and require subsequent lattice simulations [18-20] 2 
to study thermalization [26, 27] and ultimately reheating. The subsequent decay of preheat 
matter fields is seldom considered, since it requires a concrete particle physics implementation 
of inflation as in [4-7]. An observable by-product of preheating can be gravitational waves 
[28-37]. The majority of these phenomenological studies focus on one (see the reviews [1, 2]) 
or several [15, 38-42] inflaton fields, coupled to a single degree of freedom that becomes 
light at the vacuum expectation value (VEV) of the inflatons. For example, in standard, 
single-field chaotic inflation with potential V(ip) = m 2 <^ 2 /2, one can introduce another light 
scalar field that couples to the inflaton via the interaction Lagrangian C = —g 2 x 2 (p 2 /2, such 
that the effective mass of x 1S minimal at the VEV ip = 0. As the inflaton oscillates around 
the minimum, x-particles are produced in bursts whenever this point is traversed, leading to 
preheating via parametric resonance (stochastic resonance in an expanding universe). 

But why is this Extra Species Point (ESP) [43] <pesp = aligned with the minimum 
of V{ip)l There is no a priori reason why they should coincide (even though they may 
in certain models, as in [4]); thus, consequences of misalignment need to be investigated. 
Furthermore, there are several models of inflation, notably trapped inflation in one [43-45] or 
higher dimensional field spaces [46], among others [47-49] (see Sec. 2.1 for a brief overview), 
that employ a dense ESP distribution during inflation. These ESPs do not miraculously 
vanish once inflation is over, in fact, the particles associated with ESPs that are within reach 
during the oscillations of the inflatons can play the role of preheat matter fields. 

Our goal is to investigate the feasibility of preheating in the presence of one or several 
ESPs in the vicinity of the inflatons' potentials' minimum, without placing them directly at 
the VEV (in the latter case, resonances are generically suppressed, due to de-phasing of the 
inflatons [15, 38-40]). 

Focussing on more than one inflaton (M > 2) and dense ESP distributions, we find 
several distinct phases: during the first in-fall, the inflationary trajectory traverses a large 
number of ESPs at which single bursts of particle production take place. If ESPs are dense 
enough, this first infall can already be sufficient to transfer a large fraction of the energy into 
preheat matter fields before any oscillations take place. We provide analytic expressions for 
the comoving particle number, Sec. 3.3.1, which we test numerically, Sec. 3.3.2. 

The second phase covers oscillations with average amplitude that are bigger than the 
average inter ESP distance. The length and efficiency of this transitory phase is strongly 
model dependent and dominated by chance close-encounters of the trajectory with an ESP. 

In the final stage, the average amplitude of oscillations is smaller than the inter ESP 
distance, Sec. 3.4. Preheating in this phase can be very efficient if an ESP has the "right" 
distance from the origin: opposite to the earlier stages, narrow resonance is responsible for 
particle production. Based on a series of approximations, we compute analytically the loca- 
tion of resonance bands, which in turn, set the locations of ESPs where particle production 
can be operational for prolonged periods, Sec. 3.4.1; again, we compare our predictions with 
numerics, Sec. 3.4.2, finding good agreement. We also test the dependence of resonance ef- 
fects onto the mass distribution of the inflatons: while our analytic results are derived for 
simple mass ratios, they remain qualitatively valid for incommensurate masses. 

1 This perturbative decay is not complete [15] and can lead to troublesome long-lived scalar particles, as in 
the cosmological moduli problem [16, 17] 

2 Common public codes are LATTIEEASY [21], Defrost [22], PSpectRe [23], HLATTICE [24] and PyCOOL 
[25]. 
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We summarize and discuss our results in Sec. 4, paying particular attention to trapped 
inflation, since ESPs are necessarily densely distributed in this model. 

To complement the analytic and numeric results based on treating inflatons as back- 
ground fields, we incorporate back-reaction in Sec. 5 by performing lattice-simulations (based 
on a modified version of LATTICEEASY [21]) in selected cases. We find resonance to be 
robust and in line with the analytic predictions; if compared with resonances of a massive 
matter field at the VEV of the inflaton, we find resonances to be more efficient than the ap- 
proximate analytic results indicate, Sec. 5.1, high-lightening the need for lattice simulations. 

To set the stage for the above studies, we start with a brief introduction to the role of 
ESPs in multi-field inflation in Sec. 2.1. Our general setup is described in Sec. 2.2. We follow 
with a brief review of preheating with a single ESP at the VEV of the inflatons to introduce 
our notation and provide a reference point for our more complicated scenario. 

Readers familiar with ESPs may skip Sec. 2.1, while readers familiar with preheating 
may want to skip Sec. 3.1. 

2 Setup 

In this section, we motivate the presence of ESP during and after inflation and provide the 
field content used in this paper. 

2.1 The role of ESPs in Multi-Field Inflation 

An Extra Species Point (ESP) is a point, or more generally an extended location (ESL), in 
field space at which additional degrees of freedom become light; consequently, they have to 
be incorporated into the low energy effective field theory, if a trajectory in field space comes 
close to the ESP, to allow for particle production and back- reaction of the produced particles. 

Such ESPs are generic in moduli spaces of string theory, see [50-58] for some examples, 
with many applications, such as moduli trapping [43, 59-63], effects on inflation [47-49], 
trapped inflation [43-46] (see also [64-69]), modulated trapping [70, 71] and preheating, 
among others. 

These additional degrees of freedom are often associated with enhanced symmetries; for 
example, if two D-branes come close to each other (the distance between them is the modulus 
field), strings stretching between them become light (the new degree of freedom) leading to 
an enhanced gauge symmetry [50]. In essence, we are dealing with the string-Higgs effect 
[59, 72]. 

The positions of ESPs/ESLs are characteristic for the low energy effective field theory 
at hand, and little can be said in general about their distribution, which can range from 
dense, as in trapped inflation [43-46], to sparse, as in [70]. As a consequence, the effect 
of ESPs are often studied from a phenomenological point of view [43, 70], supplemented 
with concrete case studies (for example, monodromy inflation [45] is a realization of trapped 
inflation [43, 44]). 

In this paper, we model the presence of ESPs during preheating by the inclusion of 
additional scalar fields Xd with a negligible bare mass that we couple to the inflatons via 
quadratic interactions; the vector a parametrizes the position of the ESP in an AA-dimensional 
field space, see Sec. 2.2. More realistic setups could involve the inclusion of gauge fields, as 
in [73], or more general couplings, which we postpone to future studies. 

Before turning our attention to preheating, we provide a brief summary of three ex- 
emplary inflationary setups that rely on the presence of ESPs: trapped inflation [43-46], 
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grazing ESP encounters [70] and modulated trapping [70, 71] (this selection is not complete 
and chosen for illustrative purposes). In these examples, the inclusion of several ESPs (not 
aligned with the VEV of the inflatons) during preheating is not an option, but mandatory, 
motivating the distributions we consider in Sec. 2.2. However, the reader should be aware 
that ESPs are a common occurrence in almost all effective field theories and their presence 
is at the core of any preheating study - it is just assumed that their location coincides with 
the inflaton's VEV, an assumption that commonly lacks justification. 

2.1.1 Trapped Inflation 

In trapped inflation [43-46], the backreaction of the produced particles onto the inflaton(s) 
is taken into account, leading to a modification of the low energy effective potential: as the 
ESP is traversed, particles are produced which lead to an attractive force towards the ESP. 
An encounter can lead to brief temporary trapping if enough particles are produced (brief, 
since the number of particles and thus backreaction redshifts with a~ 3 ), but even if particle 
production is minor, i.e. when the ESP is grazed with a reasonably large impact parameter, 
it can still lead to observational signals, see Sec. 2.1.2. 

The term trapped inflation is commonly used if several ESPs are either encountered 
head on [43-45], as in monodromy inflation [45, 74], or generically grazed [43, 46] if the 
field space is higher dimensional. The repeated encounter with ESPs and the resulting extra 
contributions to the potential can prevent the inflatons from rolling down fast, prolonging 
inflation and consequently rendering steeper potentials admissible. If the number of inflatons 
is large J\f 3> 1, and ESPs are ubiquitous, a terminal velocity in field space results that is 
independent of the potential [46] . 

As an example, consider modelling the additional degrees of freedom by massless fields 
quadratically coupled to the inflatons, g 2 {(p — <P%Sp) Xa/^> at ESPs with average inter ESP 
distance x; further, assume a quadratic potential for all inflatons with the mass set by the 
COBE normalization and start at a distance of 0.2 M p from the origin (regular slow roll 
inflation is over at this point). If g ~ 0.1, x ~ 5 x 10 -5 M p and N > 37 (see [46] for details), 
the terminal velocity v% ~ gx 2 is small enough to drive an additional sixty e-folds of inflation 
and satisfy all observational constraints. The same mechanism can work with almost all steep 
potentials, as well as higher dimensional extra species locations [46], opening up regions of 
field space for inflationary model building, previously deemed unsuitable. 

Observational consequences are varied: the leading order contribution to correlation 
functions are altered [46] while retaining the successes of slow roll inflation: trapping effects 
slow the fields down, requiring the replacement of potential slow roll parameters with the 
Hubble slow evolution parameters, such as e ~ H/H 2 , which in turn are set by the back- 
reaction of the light particles. If the trajectory turns, additional contributions arise due 
to the interaction with isocurvature modes. Further, once the trajectory moves away from 
the ESP, the formerly light particles become heavy again. Their backscattering onto the 
inflaton condensate leads to an infra-red cascade, causing additional bump like features in 
the power-spectrum and non-Gaussianities of peculiar shape [75-78]. 

After trapped inflation terminates, either because fields start to evolve fast or because 
the energy density in the additional degrees of freedom becomes comparable to the one in 
the inflatons, preheating commences in the presence of many ESPs. 

As we shall see in this paper, if ESP distributions are dense, particle production during 
the first in-fall towards the VEV can be efficient before even a single oscillation of the inflatons 
is complete; further, any remaining energy can be drained in a subsequent resonance phase. 
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This qualitatively different type of preheating is the expected scenario for trapped inflation 
in higher dimensions [46]. 

2.1.2 Grazing ESP Encounters 

If ESPs are less dense, the inflationary trajectory may still encounter one or several of them 
during the last sixty e- folds of inflation [70]. Back-reaction generically causes the trajectory 
to bend and slow down, leading either to temporary trapping or a slightly prolonged phase 
of inflation. 

Which one occurs depends on the coupling strength to the light fields and the impact 
parameter /x, that is the distance of closest approach to the ESP of the unperturbed trajectory. 
If the value of the coupling constant depends on other fields (generic in string theory), 
modulated trapping results [71], see Sec. 2.1.3; but even if the couplings are taken to be 
constant, additional observational signatures are possible: if the coupling is strong, g ~ 1, 
and the impact parameter is not too big, temporary trapping occurs: the resulting trajectory 
is chaotic [70], and as a consequence, the duration of trapping depends sensitively on the 
impact parameter, which in turn fluctuates due to isocurvature perturbations. Since the 
time being trapped is not a smooth function of [i, the resulting fluctuations in the power- 
spectrum are usually in excess of the COBE normalization; thus, most temporary trapping 
events are ruled out during the last sixty e-folds of inflation. However, if the trajectory bends 
slightly and slows down without turning around, the dependence of the number of e-folds on 
/i, and thus isocurvature perturbations, is smooth, leading to additional contributions to the 
power-spectrum and non-Gaussianities, consistent with current bounds. In addition to these 
signatures, localised and oscillatory features are also present via IR-cascading [75-78]. 

Again, ESPs do not vanish once inflation terminates, but need to be incorporated into 
the study of preheating. As we shall see in this paper, even a single ESP at the right location 
can lead to prolonged resonances, which are cut short if the location of the ESP coincides 
with the VEV of the inflatons. 

2.1.3 Modulated Trapping 

To conclude our brief overview, we would like to comment on modulated trapping [70, 71], 
which is the application of modulated preheating [79, 80] to an ESP-encounter during infla- 
tion. If the coupling constant(s) depend on other fields, their fluctuations can be converted to 
curvature fluctuations during the encounter: instead of a dependence on isocurvature fields 
via the impact parameter, the effectiveness of trapping, or more precisely the decrease in 
speed of the inflatons, is modulated by fluctuations in additional field(s) via the coupling 
constant. 

The process of converting fluctuations is similar but not identical to a grazing ESP 
encounter, ultimately leading to similar observational signatures [70]. Naturally, both effects 
are likely to coexist in realistic models of inflation. 

To summarize, ESPs are ubiquitous in moduli spaces of string theory and depend- 
ing on their concrete distribution, a wide array of observational signatures, such as non- 
Gaussianities, are possible. In the absence of a thorough understanding of their distribution 
in the region of field space that inflaton(s) traversed during the last sixty e-folds in our uni- 
verse, it is prudent to investigate their effect onto the end of inflation, particularly preheating. 
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2.2 A Generalized Setup for Preheating 

Having motivated the presence of ESPs during/after inflation, we specify the concrete phe- 
nomenological implementation used throughout. 

2.2.1 Inflaton fields 

Consider J\f inflaton fields ipi, i = 1 . . .A/", with canonical kinetic terms and quadratic poten- 
tials 

N 1 

F(^,..., W ) = ^-m^f. (2.1) 

i=i 

This potential need not be valid during inflation - we just assume that we can expand the 
inflationary potential and diagonolize the mass matrix to describe the oscillations of the 
inflatons around the minimum. In general, the masses conform to a distribution that differs 
from mi = nij for all for example, in N-flation [81] they conform to the Marcenkov-Pastur 
distribution [82], spread around to CO be ~ 

6.2 x 10 6 /V8nG ra 1.2 x 10 6 M P which is set by 
the COBE normalization [83]. Occasionally, we set 

M% = ± = 1. (2.2) 

Since we do not focus on a particular inflationary model, we consider a simpler mass distri- 
bution: after arranging the fields according to their masses such that rrii < rrij Hi < j, we 
spread the m; around 

m = l(T 6 Mp. (2.3) 
If not stated otherwise, we distribute the masses over the interval 1 < m,/m < 1.5, that is 

m i =m(l+ 2 [^~_ 1) 1) ) = Vfcm , (2.4) 

so that mass ratios are simple fractions (enabling a simple analytic description of some 
resonance effects). 

To check for a possible dependence of results on the chosen distribution, we occasionally 

use 

_ 2 M + 2i - 3 2 _~ 2 

m% = 2(M-l) m = ' ^ ' 

as in [40], so that the average square mass is 

( ft ylf ^n 2 -lf (AA + ^-3) 2_ AA 2 +AA(AA + 1)-3AA 2 (2 6) 

i=i i=i x ' 

the lowest one is = m 2 /2, and the greatest one is m? N = 3m 2 /2 (masses are incommensu- 
rate). Note that (2.4) and (2.5) are mutually exclusive definitions of the coefficients /3j and 
Pi. One can motivate such narrow mass-distributions if (2.1) is also valid during inflation 
and all fields have comparable initial energy. In this case the lightest fields dominate the 
energy density at the end of inflation, and are thus responsible for preheating, because heavy 
fields have already relaxed to the minimum of their potential during inflation [38]. 
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We also need to specify the field values at the end of inflation, that is the initial values 
for preheating. In the absence of a concrete model, we provide each of the inflatons in (2.1) 
with the same energy at the end of inflation summing up to m?ipp r /2, where we set 

^ = |^(t pr )| = 0.193 M p , (2.7) 

in line with single- field models [8]. These choices yield 



(analogously for fit) as initial values for the inflatons. For simplicity, we further set all initial 
velocities to zero 



<Pi(t pr ) = 0, (2.9) 

since the speed towards the end of inflation is still slow roll suppressed. While we use 
these initial conditions most of the time, we occasionally deviate from them for illustrative 
purposes, for example, setting 

<Pi(tpr) = Ppr = <P0 ■ (2-10) 

We would like to point out that our conclusions turn out to be largely insensitive to the 
choice of masses as well as the initial conditions. 



2.2.2 Matter Fields 

We wish to include several preheat "matter" fields that become light at different ESPs in 
the vicinity of \tp\ = 0; we deliberately avoid putting an ESP at \<p\ = 0, since preheating of 
such a field is well understood (see [8, 9, 15, 38, 40-42, 84] for a small, non-representative 
collection, with particular focus on multiple inflatons) and generally suppressed if N > few 
[15, 38-40]. 

We consider three types of ESP-distributions: 

1. Along a grid 

<esp=Q + «^ (2.11) 

with O) £ Z. This distribution is almost always used for numerics in the absence of 
back-reaction. 

2. A random distribution, placing ESPs randomly over the accessible field space up until a 
desired average inter ESP distance is reached; used to check if the grid-like distribution 
leads to artefacts (it does not). 

3. A Gaussian distribution centred at (p = with variance a and desired x near the 
origin; used for lattice simulations, focussing on the late stages of preheating when the 
oscillation amplitude is smaller than a. 
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For simplicity, we ignore interactions between matter fields and their bare masses (a 
large bare mass suppresses resonances); we further assume that they couple with the same 
strength 



g = 5 x 10" 



(2.121 



to the inflatons. This coupling is chosen to enable easy comparison with prior preheating 
studies such as [8, 15, 40]. Note that this coupling is often assumed to be much bigger for 
models that employ ESPs during inflation, in order to yield measurable signatures [43, 44, 46, 
70, 75-78]. The qualitative features of preheating discussed in this article are not contingent 
on a finely tuned g and remain valid in more general cases. 
All in all, the complete Lagrangian reads 



2 2 



1=1 



i=i 



v?. 



(2.13) 



To discuss preheating, it is useful to define a non-adiabaticity parameter for a given ESP 



LOT 



g 2 0(0 - <^ ES p) - Hk 2 /a 2 



(g 2 (0- <a=sp) 2 + k 2 /a 



2\3/2 



(2.14) 
(2.15) 



In most figures, time is given in terms of the number of oscillations in a hypothetical field 
with mass m = 10 _6 Mp, that is 



N 



mt 
2^' 



(2.16) 



which should not be confused with the number of fields J\f or the number of e- folds (not used 
in this paper). 



3 Preheating in the Absence of Back-reaction 

Treating the evolution of the inflaton fields as a background for the matter fields, that is 
ignoring back-reaction, we would like to investigate the feasibility of preheating in the setup 
of Sec. 2, with the field content and interactions of Sec. 2.2. Concretely, we would like to 
examine the effectiveness of particle production at ESPs during three distinct phases: 

1. The First In-fall: directly after inflation, the inflaton fields evolve rapidly towards 
0=0, traversing many ESPs along the way. 

2. Large Amplitude Oscillations around 0=0 follow the first in-fall, with A 3> x. 

3. A Resonance Phase with small amplitude oscillations around 0=0, with A -C x. 
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For the important first and third phase, an analytic description of preheating is possible; 
the brief intermediate phase is strongly model dependent with at most intermittent particle 
production. Therefore, we focus on the first and third phase. 

Analytic results are compared with numerics, entailing the integration of the equations 
of motion 

dV 

<Pi + 3H<pi + — = 0, (3.1) 
dipt 

xf + 3Hxt + u;lxi = 0, (3.2) 
with the time dependent oscillation frequency 

ul = ^ + g 2 {0-rE S p) 2 i (3-3) 
in conjunction with the Friedmann equation 

V i 

using common C++-routines. These results are supplemented by some computationally 
more demanding lattice simulations in Sec. 5, which incorporate back-reaction of the preheat 
matter fields onto the inflatons. 

To set our notation and provide a self contained treatment, we review some well known 
properties of preheating at a single ESP in Sec. 3.1, following primarily [8, 38, 85]. 

3.1 Review of Preheating at a Single ESP 
3.1.1 Inflaton Dynamics during Preheating 

Introducing the angle variables Qi, defined via 



0i ~ V 4^ H C ° S 9i ' = V 4v^ F &i ' (3 ' 5) 

the background equations of motion become 

H = -3# 2 -^ Yj cos2 9i ' ( 3 - 6 ) 



3 

$i = rrii + -H sin 20*. (3.7) 

During most of the preheating phase m; > H holds, so that Oi ~ nut up to some irrelevant 
phase. Thus H ~ 2/(3t) to leading order, so that a oc t 2 / 3 and 

(fi(t) sin rriit = & smrriit , (3-8) 

VS-irAfrriit 

that is, the fields oscillate with ever decreasing amplitude <pi{t) oc 1/t. Further, even if fields 
start out in phase, they quickly de-phase unless the masses are identical. 
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Figure 1. The first two instability bands (dark blue) of the Mathieu-equation for small values of q 
(computed analytically for q <C 1). A k < is not possible during preheating. 

In the following, we need the average distance of the trajectory from the origin (p, as 
well as the average velocity v (averaged over one oscillation). With the above approximation 
we get 



where we used the virial theorem to estimate the velocity. 

3.1.2 Parametric Resonance at a Single ESP Driven by a Single Inflaton 

Consider a single ESP at tp = e so that Ci n t = —g 2 (f — e) 2 x 2 /2. Since the oscillation 
frequency is time dependent, resonances can take place, increasing the particle number 

Pk 1 _ w fc (\Xk\ 2 . | ,2\ 1 



n k = — - 1 = — ^ + \Xkr)- 1 (3.11) 

explosively. To encounter resonance, consider <p> <C e, ignore the expansion of the universe 
and define = 4(/c 2 + g 2 e 2 )/m 2 , q = Ag 2 ecj)/m 2 and 2z = mt + 7r/2 so that (3.2) becomes 
the Mathieu-equation 

xl + {A k -2qco S 2z)xk = 0, (3.12) 

where a prime denotes a derivative with respect to z. According to the Floquet-theorem, 
solutions to (3.12) are of the form 

Xk = rj{z) exp^z) , (3.13) 

where rj is periodic with period tt and the Floquet index (i^ (in more general cases, we simply 
refer to p k as a characteristic exponent) is positive in certain bands, see Fig. 1; for example, 
in the first band for q <C 1 one finds [86] 
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1 2 3 4 5 6 7 8 0.5 1 1.5 2 2.5 

t [N] t [N] 

Figure 2. Numerical solution of the particle number for narrow (left, (Aj.;q) = (1;0.2)) and broad 
(right, (Ak\q) = (400; 200)) resonance, ignoring the expansion of the universe. Time is plotted in 
units of N = mt/(2n) with m = l0~ 6 M p . 



As a consequence, resonances occur for m(l — q/2)/2 < k < m(l + q/2)/2. For 99 > g 2 e/(32ir) 
this narrow resonance dominates over classical reheating effects [8]. 

In an expanding universe, modes do not stay within a given instability band, but traverse 
through them to approach the stable region near A k = 0, q <C 1. As a consequence, the 
effectiveness of preheating is reduced: for example, after At ~ qH -1 a mode in the first 
instability band will have left said band, leading to a net amplification of the particle number 
by ~ enp(q 2 m/(2H)). Thus, a condition for effective preheating is 

q 2 m>2H, (3.15) 

that is ip > my / mH/2/(2g 2 e). 

The above description is valid towards the end of preheating, when the oscillation am- 
plitude is small. At earlier times eft < Mp and we need to consider (f> ^> e, so that 

Xk + {k 2 +g 2 ip 2 ) X k = 0, (3.16) 

if expansion effects are ignored. Defining A^ = 2q + k 2 /m 2 , q = g 2 (p 2 / (4m 2 ) and z = mt, we 
arrive again at the Mathieu equation, leading to broad resonance, see Fig. 2. 

Particle production in the broad resonance regime takes place in bursts, whenever ip ~ 
so that the effective mass of the matter field is light. To be concrete, if the non-adiabaticity 
parameter in (2.14) becomes large 

T)>1, (3.17) 

particle production can take place. 

To include the expansion of the universe, it is convenient to define X = a?l 2 Xk s ° that 
the equation of motion becomes 

X k + i4X k = 0, (3.18) 

with u 2 . ~ k 2 /a? + g 2 (f> 2 sin 2 (mt) and a cx i 2 / 3 . The resulting solutions still show ample 
particle production that appears to be stochastic, see Fig. 3, and carries therefore the name 
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Figure 3. Numerical solution of the particle number for the k = Am mode during stochastic resonance, 
that is broad resonance in an expanding universe, with g = 5- 1Q- 4 and m = M P . 



stochastic resonance. It should be noted that this process is fully deterministic, and can be 
understood analytically [8]: in the WKB approximation, the adiabatic solutions (rj <C 1) to 
(3.18) are 

X k (t) = -^ =e i J" w * di + J^ =e - < J" w * dl , (3.19) 

where the time dependent Bogolioubov coefficients have to satisfy |a/c| 2 — \/3k\ 2 = 1- The 
comoving particle number is n& = |/3fc| 2 . Expanding <p(t) around the j's particle production 
event and defining r = fc*(i — tj) as well as k = kj (ak+) yields 

^ + ( K 2 + r 2 )X fc = 0, (3.20) 

which describes scattering at parabolic potentials. A computation of the transfer-coefficients 
relates the Bogolioubov coefficients before and after the particle production event (while 
•q S> 1), leading to [8] 



e-™ 2 + (1 + 2e"™ 2 K " 2e-™ 2 / 2 V / l + e-™V n i( 1 + 4) sin C > ( 3 - 21 



where Q\ ot is a time and fc-dependent total phase. Thus, the step like change of the particle 
number as well as the occasional drop can be understood analytically (a drop takes place if 



sm0L > e~< I 2 ) 



- 12 - 



50 



40 



+ 30 



S 



20 



10 




Evolution of the k — mode for different numbers of infiatons J\f with e — and 
Resonaces are suppressed as M is increased due to the dephasing of fields (inflaton 
masses are chosen according to (2.6)). 



Figure 4. 

a = 5- 10" 



As evident, resonances can efficiently increase the particle number in an expanding 
universe. Naturally, once the particle number is large and the energy density in the matter 
field becomes comparable to the energy of the inflaton, one needs to incorporate back-reaction 
effects. However, to decide whether or not parametric resonance is important, the simplified 
treatment of this section is fully satisfactory. 

3.1.3 Parametric Resonance at a Single ESP Driven by Several Infiatons 

How does preheating change if several infiatons are present? A straightforward generalization 
of the interaction Lagrangian to 



-int 



2, ,2 



(3.22) 



yields an, at first sight, surprising result for |e| = 0: resonances are strongly suppressed if 
more than a few infiatons with different masses are present [15, 38-40], see Fig. 4. However, 
this can be understood easily: infiatons generically run out of phase due to differences in 
their masses. As a consequence, the trajectory never crosses the position where the ESP is 
located (the origin), the non-adiabaticity parameter remains small and particle production 
is suppressed. 

In case of two infiatons, an additional effect can keep the efficiency up to some degree: 
if the mass ratio of the infiatons is not a rational number, the driving force is only quasi 
periodic. As a consequence, stability bands dissolve into a no- where dense set, leading to 
a positive generalized Floquet index for almost all modes [41, 42] (Cantor preheating); this 
result is based on spectral theory [87-93], and numerics. Unfortunately, the amplitude of the 
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Figure 5. Schematic of an ESP-encounter in a two dimensional field space. If the trajectory comes 
close to the ESP, fx < \i c = {\0\/g) 1 ^ 2 , the adiabaticity condition r\ < 1 can be violated, and particle 
production can take place. 



characteristic exponents can not be computed analytically and it appears that most of them 
are generically small; as a consequence, the suppression due to de-phasing wins if more than 
two inflatons are included [38, 40]. 

3.2 Overview: The Three Phases of Preheating 

In the last section, we discussed preheating in the presence of a single ESP at the origin: to 
allow for explosive particle production, we saw in (3.17) that 77 > 1 is a sufficient condition 
(77 is the non-adiabaticity parameter in (2.14) and we considered small /c-modes). During 
encounters with a non-zero impact parameter /x (see Fig. 5) this condition becomes 

H<Hc (3.23) 
where we defined the critical impact parameter 




Hc=\—. (3-24) 



Even though v = \<p\ is slow roll suppressed during inflation, particle production may still be 
strong enough to lead to trapped inflation [43-46] or additional contributions to observables 
[70], if ESPs are dense or an encounter is close (see Sec. 2.1.1 and 2.1.2). During preheating, 
the critical impact parameter increases sharply, making ESPs in wider regions of field space 
relevant for particle production. 

If we expand the trajectory near the ESP encounter (p ~ </? E sp + ft — vt, we can write 
the equations of motion for as 

^ + ^ + r 2 )X k = 0, (3.25) 
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Figure 6. Post inflationary trajectory in a two-dimensional field space, (p%/M p over ipi/M„ with 
initial conditions <fii(to) = ^2(^0) = 0.2 M p , ipi(to) = (^2(^0) = and masses mi = 10~ 6 Mp and 
ni2 = 1.5 • lCP 6 Afp. The dotted blue circles denote regions that the trajectory does not leave on 
average after t/N — 0.5,8,16, corresponding to the three phases: first in-fall (t/N < 0.5), large 
amplitude oscillations (0.5 < t/N < 16) and resonance (16 < t/N < 50). The corresponding total 
particle number in the k — 0-mode of a grid-ESP distribution with inter ESP distance x = 0.0033 M p 
is plotted in Fig. 7. 



with r = ^fgvt and 



gv \a A 1 



(3.25) is formally identical to (3.20), with the above redefinitions of r and k. Therefore all 
prior results, such as the particle number n 3 k in (3.21), carry over (the index j now counts 
close encounters with the same ESP, as opposed to oscillations of the inflaton). 

Equipped with this knowledge, we are ready to discuss particle production for actual 
trajectories: in Fig. 6 we plot an exemplary trajectory in a two dimensional field space, 
with the corresponding particle production at ESPs distributed along a grid with inter ESP 
distance x = 0.003 M p in Fig. 7. 

We initialize (/'l^o) = ^2(^0) = 0.2 M p , ^i(^o) = ^2(^0) = 0, and follow the fields 
evolution in an expanding universe. We chose zero initial speed to mimic the low speed 
at the end of inflation. In the absence of large Hubble friction, the speed grows quickly, 
increasing the critical impact parameter. Since /U c 3> x, many ESPs are within reach during 
the first in-fall and a seizable quantity of particles can be produced (we lump all particle 
species together in the particle number depicted in Fig. 7), accounting for the increase to 
ln(nfc + 1) ~ 7 in the insert 3 . Even though the particle number at each individual ESP is 
low because each one is only passed once, the resulting overall increase in n^, and thus the 
corresponding energy transfer, can be considerable. The end of this phase, which is absent in 
the single ESP case, is marked by the blue dotted circle on the left of Fig. 6 and the vertical 
line in the insert of Fig. 7. 



3 If ESPs are dense and M is large (larger than in this section), back-reaction can lead to trapped inflation 
at a terminal velocity even for the sub-Planckian field values considered here (see the worked example in Sec. 
III.C of [46]). Once inflation ends, a seizable amount of x-particles are already present in this scenario; the 
subsequent evolution is identical to the one discussed above. 
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Figure 7. Occupation number of the k = mode in an ESP distribution along a grid with inter-ESP 
distance x ~ 0.0033 M p with two inflatons evolving as in Fig. 6. Solid vertical lines separate the 
phases of preheating, as in Fig. 6. 



During the subsequent large amplitude oscillations the fields slow down v oc 1/t, fi c 
decreases and particle production takes place only sporadically in bursts, reminiscent of 
stochastic resonance: in our example, rik stays roughly constant until t/N = 8; bursts become 
more frequent thereafter, leading up to the resonance phase. The duration and effectiveness 
of this intermediate phase is model dependent and somewhat unpredictable; it ends once the 
amplitude of oscillations falls below the inter-ESP distance 4> < x, indicated by the small 
blue circle on the right side of Fig. 6 and the vertical line at t/N = 16 in Fig. 7. 

During the third resonance phase, the oscillation amplitude of the inflatons is small. 
This phase is reminiscent of the narrow resonance regime, inducing a steady increase of the 
particle number; as we show in Sec. 3.4, particles are produced at a few ESPs that have 
the "right" location to enable prolonged resonances. Particle production ends, once the 
amplitude/speed become too small, around t/N = 50 in our example (see Sec. 4.2), and the 
particle number remains constant thereafter. 

The three phases, especially the first in-fall and the resonance phase, are generic to 
multi- field preheating in the presence of dense ESP distributions, x <C M p . We provide the 
analytic explanation of these phases in Sec. 3.3 and Sec. 3.4, which contain the main new 
insights gained in this paper (see Sec. 4.1 and 4.2 for a summary of the results). 

Efficient preheating requires that the energy density in the different matter particle 
species becomes comparable to the one in the inflatons before the resonance phase ends. In 
this case one has to include back-reaction, see Sec. 5, as well as the decay of preheat matter 
fields [4, 5] (not covered in this paper). 

An important result of such multi-field preheating is the presence of two distinct particle 
classes: those produced at many ESPs with low occupation number (first in-fall and inter- 
mediate phase) and particles produced in great number at a few ESPs (resonance phase); it 
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Figure 8. Schematic of a trajectory during the first in-fall. ESPs in the green shaded regions are close 
enough to allow for the violation of the adiabaticity condition r\ < 1, and thus particle production. 
(J>c = {W\l g) 1 ^ 2 is the impact parameter up to which we consider ESPs (/i rs /i c for most ESPs in the 
green region; this approximation becomes better, the larger Af is). 



is tempting to identify the subsequent decay products of particles in these two classes with 
dark matter and standard model particles respectively. However, one would need to consider 
a specific scenario to concretize this speculation, which goes beyond the scope of this study. 

3.3 The First In-fall 
3.3.1 Analytic Description 

During the first in-fall, many ESPs are close enough to the trajectory to allow for adiabaticity 
to be violated and particles to be produced, see Fig. 8 for a schematic. We would like to 
compute the occupation number in the k = mode for comparison with numerics, as well as 
the total particle number 

n x = J n k d 3 k, (3.27) 

where we combine all the different matter species. We closely follow [46] in this section. 
Since x-nelds start out in the vacuum, only the first term in (3.21) contributes to the particle 
numbers during the first in-fall, 




(3.28) 



where fj,$ is the impact parameter of the ESP at <p^ SP , we used (3.26) and ignored expansion 
effects (a = 1), since inflation already terminated. To compute the total particle number, we 
replace the discrete ESP distribution by a homogeneous ESP density p ESP = x~-^ so that 

dn fc = {n k ) p ESP dV (3.29) 
= (n fe )p ES pVV-id^, (3.30) 

where V/v-i is the volume of an M — 1 dimensional hyper-sphere, and (nk) denotes the 
average particle number at ESPs in the infinitesmal volume element dV shown in Fig. 8. As 
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Figure 9. Left: comparison of (no) with the large TV-limit exp(— it) (taking /j, = fi c for all ESPs in 
dV). Right: velocity-integral J- (AT, a) for different values of a with the corresponding fits. 



a rough estimate for the radius of the hyper-sphere, we take the critical impact parameter 
Mc = y/v/g, since particle production at ESPs further away is strongly suppressed 4 , leading 
to 



dn 



r( 



(n k ) dp , 



(3.31) 



where T is the Gamma function. 

Let us focus on the k = mode to see how rapidly (no) approaches the limit value 
exp(— 7r) as TV is increased. Based on the definition of the average, we get 



/exp 



Tfl ,,2 



TV-1 



IT 



Af-1 

-— r 



TV-l 



, 7T 



(3.32) 



(3.33) 



which we plot in the left panel of Fig. 9. Evidently, the limit is approached rapidly, providing 
an excellent approximation for TV as low as 20, while still giving the right order of magnitude 
for lower TV. This is not surprising, since almost all of the volume in a sphere is near the 
surface if TV ^ 1, so that the impact parameter to almost all relevant ESPs is close to u c . 

Using this approximation, we can compute the total occupation number in the zero- 
mode by integration over the first in-fall to 

ft* 



71, 



tot 



dn = Pesp (no) • / V/V-i v dt 



M-i r 

2 



r(^i) 



where we defined the velocity integral 

T{Af, a) 



V 2 " 



"(t) dt. 



(3.34) 
(3.35) 

(3.36) 



o 



4 See [46] for an improved estimate in the large Af limit; the simpler estimate above is satisfactory for our 
purposes (we slightly underestimate particle production). 
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a 


Cl 


C2 


0.5 


-7.6777(4) 


6.360(5) 


2 


-7.653(4) 


17.95(4) 



Table 1. Coefficients of the approximation in (3.37) for the velocity-integral in (3.36) via linear 
regression. 



The latter can usually be approximated by 

F(N,a)^e Cl{a)+(yAa)N , (3.37) 

with constant coefficients ci, ci- For example, using the mass distribution of Eq. 2.4 and the 
initial conditions specified in Sec. 3.2 for all of the M inflatons leads to the coefficients in 
Tab. 1. The corresponding regression is shown in the right panel of Fig. 9. 

The analytic expression of n ot is easily compared to numerical results, providing a check 
for the validity of the approximations, see Sec. 3.3.2. 

It is instructive to plot the total co-moving particle number in the zero mode at all 
ESPs after the first in-fall over the number of fields N for varying inter-ESP distances x, see 
Fig. 10: the solid black line corresponds to the particle number no generated at the end of the 
entire preheating phase with M = 1 and a single ESP at the origin (the canonical example 
of preheating); thus, this line serves as an indicator whether or not preheating is efficient. 
As evident, for dense ESP distributions and many inflatons, x < 0.0057M p and N > 40, 
more particles are produced during the first in-fall than during the entire preheating phase 
in the corresponding single field model. If ESPs are denser, less inflatons suffice. Thus, the 
first in-fall alone has the potential to compensate for the (potentially) decreased efficiency 
of preheating due to de-phasing effects [15, 38, 40]. Furthermore, as we shall see in Sec. 3.4, 
the late resonance phase can be prolonged if an ESP happens to be at the "right" distance 
from the origin. 

Since the total comoving particle number n x is related to the energy in the x- sec ^ OT i 
and thus better suited to gauge the effectiveness of preheating, we would also like to provide 
an analytic approximation for n x . Using the definition of n k in (3.28) we get 

(n k ) = (n ) e , (3.38) 



so that 



/ tot *^ f 7 2 tot J 7 

' h ~- 1 Uk Wf = Wf Jo Hk 

^9 — ^ /jy-n x I v ^ \ I k e 9 dk ) dt - ( 3 ' 39 ) 



(2tt)3 3 r (^=i) 



vo 



With the definition of the velocity-integral, we arrive at 



F(^) ' (2vr)3 



% = 9~ ~r x~ N ■ . (3.40) 
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Figure 10. The total comoving occupation number after the first in-flall over the number of fields 
JV for different inter ESP distances x. The solid black line corresponds to no a t the end of preheating 
in presence of a single field and an ESP at the origin, see Fig. 4. 

(3.40) provides an analytic understanding of particle production during the first in-fall and 
is our first main result. In its derivation, we made a series of approximations: we ignored 
expansion effects (justified since the first in-fall typically takes much less than a Hubble 
time), did not incorporate back-reaction (valid as long as occupation numbers do not grow 
too much) and considered particle production at ESPs up to the critical impact parameter 
Mc = y/g/v only. 

Naturally, to gauge the ultimate effectiveness of preheating, back-reaction effects need 
to be taken into account once particle numbers in the x~ ne lds are high (lattice simulations 
are needed) or if the corrections to the effective potential brought forth by the produced 
particles rivals that of the bare one 5 (the WKB approximation can be used). 

3.3.2 Comparison with Numerics 

We would like to test the validity of the analytic expressions in (3.35) and (3.40), which 
required a series of approximations. 

To this end, we compute numerically the comoving occupation number in the zero mode 
riQ * at the end of the first in-fall, t = t*, for different inter-ESP distances x and distributions 
(grid like and random, see Sec. 2.2.2) for up to four inflatons N = 2, 3, 4 (for N = 4, we only 
consider inter ESP distances down to x = 0.02 M p due to increasing computational demand), 
see Fig. 11. 

For example, if inflation is driven at the terminal velocity in a higher dimensional field space, as in [46], 
back-reaction is important throughout; further, inflation ends once p x becomes comparable to the energy in 
the inflatons, leading to a smooth transition to preheating. 
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Figure 11. Comoving particle number in the k — mode, ntf*, after the first in-flall for different 
inter-ESP distances x and number of inflatons M . Solid lines: analytic approximation in (3.35); filled 
circles: numerical result with ESPs distributed along a rectangular lattice, see Sec. 2.2.2; open circles: 
random ESP distribution with the same average inter ESP distance as the grid. 

We observe that the analytic approximation does a good job in explaining particle 
production during the first in-fall. The scatter around the analytic result is caused primarily 
by evaluating n^ 1 at a fixed time, which is sensitive to the small spikes in Fig. 7 (one could 
average over the first plateau to reduce the scatter). These spikes are caused by chance, close 
ESP encounters. A concrete example is present for a regular grid with x = 0.016 in Fig. 11: 
our initial conditions put the trajectory right on top of an ESP, leading to increased particle 
production (as a consequence, the corresponding data-point falls out of the range depicted 
in Fig. 11). Of course, such chance encounters can not be recovered by our analytic result. 

To check whether or not ESP positioning along a grid influences the outcome quali- 
tatively, we also performed numerical runs with randomly placed ESPs that have the same 
average inter-ESP distance in the region of field space of interest. Discrepancies are small 
and of the same order as the scatter around the analytic expression, which is expected. 

Changes in the initial conditions lead to altered analytical as well as numerical results, 
for two reasons: firstly, choosing a different starting point can prolong/shorten the first in- 
fall. Secondly, changing the initial speed alters the value of the velocity integrals T: an 
increased speed leads to a larger fi c and thus boosts particle production. Both effects are 
evident in Fig. 12, but they are ultimately minor, leaving the qualitative features unchanged. 

Thus, the analytic expressions in Sec. 3.3 provide an excellent description of particle 
production during the first in-fall as long as back-reaction effects are negligible. 
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Figure 12. Dependence of particle production during the first in- fall on initial conditions. Left: 
comoving occupation number is plotted over time for Af = 2, x = 0, 005 Mp and several initial 
conditions (ipo; ipo) in units of (Mp; 10~ 6 AfJ>) for all inflatons, except in the third and forth case 
from above, for which we also set ^i(O) = 0, 2 Mp and ip%(0) = 0. Right: The corresponding value of 
the velocity integral J-{N , 0.5), compared to the one resulting from our "standard" initial conditions. 
Changes in particle production are present (as expected) but minor. 



3.4 The Resonance Phase 

After a brief intermediate phase, the oscillation amplitude is smaller than the inter-ESP 
distance x, and resonances can occur. Opposite to the first in-fall, these resonances take 
place at a few ESPs only. Thus, let us focus at an exemplary ESP located at 

<Pesp = ^, (3.41) 

and investigate how much particle production in the k = mode occurs until preheating 
terminates, i.e. when the final plateau in Fig. 6 is reached. Note that the distance to the 
origin is 

d=VJf^. (3.42) 

To get an overview of possible outcomes, we first compute uq numerically, as in the 
previous section, for different values of x, see Fig. 13. The plotted value is the result of 
performing a linear fit over the final plateau (the interval is chosen to be long enough so 
that the average slope is close to zero). The top plots are for the mass distribution in (2.4) 
(simple mass ratios), while the bottom one follows (2.5) (incommensurate masses). 

Our intuition based on Sec. 3.3 indicates that an increasing x should result in decreased 
particle production: the impact parameter to the ESP is on average larger, resulting in 
a suppressed particle number according to (3.28). This intuition is wrong, as evident in 
Fig. 13: increasing x from zero leads to a steady increase of particle production with spikes 
for certain values of x. These resonances appear to be more frequent and closer to x = 
as the dimensionality of field space is increased, but also weaker, in line with the results of 
Sec. 3.1.3. Only for large values of a; is a steady suppression evident, in accord with ones 
intuition. 
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Figure 13. Comoving occupation number of the k = mode over the distance of a single ESP 
from the origin, (3.41), for two and three inflatons. Top: masses are distributed according to (2.4), 
such that ratios are simple fractions, i.e. for Af = 2 we have 1112/ mi =3/2 and for Af = 3 we get 
m 2 /mi — 5/4, m 3 /mi = 3/2. Bottom: masses are distributed according to (2.5), such that ratios are 
irrational, i.e. for Af = 2 we have m^/mx = v3 and for Af = 3 we get m^/mi = y2, m^/mi = v3- 
The first peak at x = 0.00141 M p in the top-left panel is truncated. 



Our intuition is flawed, because it is based on adiabaticity and violation thereof, i.e. we 
argued that 77 <C 1 needs to be violated for particle production to be efficient. This point 
of view is valid if the effective mass of the matter field = g\(p — v?esp|> an d thus the 

oscillation frequency u>k, is considerably larger than the inflaton masses; this is only the case 

in the early stages of preheating: the preheat matter field oscillates many times during a single 

1/2 

oscillation of even the heaviest inflaton field; since \xk\ oc uj^J , see (3.19), the particle number 
can only change during the brief moments when \<p— (p ES p\ < Mc- However, once the amplitude 
of the effective mass becomes comparable to the inflaton masses, it is impossible to separate 
the evolution of into adiabatic/non-adiabatic regimes (exponential increases in particle 
number are possible during the latter, as seen in Sec. 3.1.2). Since stochastic resonance arises 
from the interplay of adiabatic and non-adiabatic regimes, it can't be operational here and 
we conclude that narrow resonance must be responsible for the spikes visible in Fig. 13. 

We confirm this expectation in the next subsections, and, in the case of simple mass ra- 
tios, compute the location of these spikes as well as the generalized Floquet index responsible 
for the particle number's exponential increase. 

This description remains qualitatively valid for incommensurable masses, but we will 
not be able to provide exact analytic expressions. 
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3.4.1 Instability Bands 

To perform a stability analysis for X k in case of simple mass ratios, (2.4), we ignore the 
expansion of the universe so that 

N 

k 2 + g 2 Yl M ) cos ( m ^) - - 



X'k H n 



i=l 



X k = 0, (3.43) 



where we introduced z = m\t, a prime denotes a derivative with respect to z, we used 
(fi(t) = 4>q cos(mjt) (static universe) and put the ESP at <pesp = x/2. Defining 

^ = 4mf ' «1 ="^f ^ =^2-> ^ 

we arrive at an equation similar to the Mathhieu-equation 

X'; + u 2 k X k = Q (3.45) 

with 



A k - 2 £ cos ( 2 ^z) - 2 £ g« cos ; (3.46) 



The differences to the Mathieu-equation are the presence of several oscillators and the 
appearance of new parameters, and ■ As a consequence, we can not directly apply 
Floquet-theory even in the narrow resonance regime. Nevertheless, it is illuminating to 
investigate the instability regions dependence on (f>^ to decide whether or not resonance are 

(i) (i) —1 

likely. Since 4>o are decreasing in an expanding universe, <pQ oc t , a certain mode will not 
remain in an instability region, but scan several bands - knowledge of these bands enables 
us to predict the occurrence and length of preheating [94]. 

The main difference to prior studies such as [8] is the position of the ESP (it is not at the 
origin any more), leading to qualitative differences. Its main effect is the addition of a new 
term oc xj in A k (indistinguishable from a bare x _mass ) > as we h as the appearance of q2 oc xi 
(unlike a bare mass). Thus, conclusions in the small qi limit carry over directly to multi- field 
preheating with a massive x-held that becomes light, but not massless, at the origin (see 
Sec. 5.1). Surprisingly, this case has not received much attention in the literature either (see 
however [94] ) , primarily based on the intuition that a bare mass suppresses broad resonances 
(it does) and that narrow resonances are always irrelevant in an expanding universe (they 
are not, as we shall see below). 

Let us examine the case qf , q^ *C A k , that is small oscillation amplitudes and an ESP 
close to (but not at) the origin. Note that even though the presence of the ESP breaks 
spherical symmetry, it is regained in this regime since 



Ak 2 + g 2 AA(2^g + x 2 ) 
4m| 

depends only on the distance of the ESP to the origin, 



- = \(pESP\ 



E(|) 2 . (3-«) 
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and the total oscillation amplitude 



,(i)2 



V E^ 2 - (3-49) 
\ i=i 



Since we work in the limit of small amplitudes and we are particularly interested in the case 
<po < x/y/2 and low k, we approximate 



Ak 2 + g 2 Mx 2 



A k « (3.50) 



to remain consistent. We comment on the breaking of spherical symmetry by the presence 
of in Sec. 5.2. 

Defining s\ = and S2 = X' k , we get a simple system of first order differential equations 
^-s = M{z)s with M(z)=( 2/ N n ) • (3-51) 

Given initial conditions sq, its solution s(z) = U(z, zq)s(zq) is given by a linear evolution 
operator U(z,zo) that solves 

^, Z o) = MWW(2,^), (3.52) 

and satisfies U(zz,Zo) = U(z2, z%) U{z\, Zq). 

Let us denote the smallest common multiple of Tj = 27rmi / mj , that is the period of uj 2 . , 
by T. For the mass distribution in (2.4) we get 6 

T = 4tt(AA-1). (3.53) 

If mass ratios are incommensurable, oof. is only quasi-periodic and this period does not exist. 
However, the presence of spikes in Fig. 13 is not contingent on T being the exact period, 
but merely the time-scale after which oscillators tend to be (approximately) in phase again; 
practically, it appears that truncating the mass ratio to one digit and computing the resulting 
period provides a good estimate of this time-scale for our purposes 7 . 

Given this "period" and keeping in mind that the following is only exact for commen- 
surable masses as in (2.4), we get for j E N 

s(jT + zq) = (U T (zo)ys(z ) , (3.54) 

where we denned Ut(zo) =U(zq + T,zq)\ based on this property, we can perform a simple 
stability analysis, focussing on the eigenvalues of IA: if they are imaginary, solutions are 
stable, while real eigenvalues indicate an instability. The eigenvalues are 

A± = T,U T ± V(Tr^ T ) 2 -4det^ T (g g5) 



6 Only the frequencies of need to be considered, since they are (up to a factor of 2) identical to the ones 
of qi. 

7 A precise definition of T for incommensurable masses is of course desirable, but we were not able to find 
one. In all cases that we investigated numerically, such a relevant time-scale T clearly exists. 



- 25 - 



According to Jacobi's identity for the derivative of a determinant, we conclude that det Ut = 
const., since the trace of Ai(z) vanishes. Further, since Ut{z$,zq) is the identity matrix, we 
find det Ut = 1, from which we can deduce the stability condition 

\TiK T \<2. (3.56) 

In general, it is impossible to determine the evolution operator analytically; however, for 
negligible small the equations of motion are the ones of a harmonic oscillator with constant 
frequency \fA~k~ and known solutions. Since the trace is invariant under a change of basis, 
we can choose the linearly independent initial value problems s(0) = (1, 0) and s(0) = (0, 1), 
leading to 

/ cos (\/Au z) —^r= sin ( \TKi z) \ 
U(z,0) = ^ V k > yL \ M . (3.57) 

\-yJA k sin [y/A k z) cos(V^fc^) J 

According to (3.56), the transition from stability to instability regions occurs at 



cos [J AhT 



mr 



2 



1 O A k = [—) , n G N . (3.58) 

(i) (i) 

The limit of negligible q\ ,q\ becomes a good approximation in the late stages of preheating, 
since they are redshifted by the expansion of the universe. Thus, the asymptotic value of 
A k , which is set primarily by the ESPs position parametrized by x (momenta also redshift), 
determines the effectiveness of preheating and we expect the strongest particle production 

at 8 

52A k=f^ 2 ^ = ^%~- (3.59) 



4m? n IT/ " g^f T 

As we argued above, the above is valid for arbitrary positions of the ESP at distance d n — 
VJ7 x n /2 from the origin, since only this distance entered in our derivation. 

For the concrete mass distribution in (2.4), that is for masses distributed in the intervall 
[1, 1.5] • 1(T 6 M P , we get for Af = 2, 3 

x n (M = 2) =7,071 -10" 4 Mp-n, (3.60) 
x n (N = 3) = 2, 887 • 10~ 4 Mp-n. (3.61) 

In the interval 0.001M p to 0.002Af p we find indeed peaks at these positions, see top panels in 
Fig. 13, but not for smaller or larger values. This is not particularly surprising, since we made 
some harsh approximations along the way. For incommensurate masses, we find qualitatively 
the same behaviour, given that T is not identified with an exact period, but a characteristic 
time-scale common to the oscillators, as explained above. Unfortunately, the amplitude of 
the peaks is not as easily computed analytically, requiring a numerical integration. 

3.4.2 Numerical Results 

To go beyond analytic estimates, we compute numerically the characteristic exponent /x of the 
zero mode, k = 0, for varying distances of the ESP from the origin, x, oscillation amplitude of 



The index n labels the n'th distance of the ESP from the origin at which particle production is efficient, 
d n — \/J7x n /2; the x n should not be confused with the coordinates of the ESPs position in field space, Xi 
where i — 1 . . . A/". 
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Figure 14. The characteristic exponent \i of the k = 0-modc for different ESP distances from the 
origin, plotted over the initial values of the inflatons (fto (top) and time (bottom). 



the inflatons, (fto < x, and mass ratios in a static universe. We chose (ft ■ = (fto and Xi = x for 
all % = 1 . . . A/", keeping in mind that results are valid more generally since spherical symmetry 
is recovered in the small amplitude limit. 

To compute \x we perform a linear fit to ln(nfc + 1) whose slope, divided by 2ir, provides 
the desired exponent. The appearance of resonance bands (positive fi) in plots of fj, over (fto, 
see Fig. 14 to Fig. 16, indicates that the zero mode gets amplified; in an expanding universe, 
the amplitude (fto decreases over time, such that a mode is amplified for a limited time only 
while the oscillation amplitude lies within an instability band. 

In Fig. 14 we increase the distance of the ESP from the origin by changing x from x = 
over x = 0.00015 M p to 0.000515 M p . We observe the appearance of more instability bands 
as the ESP is shifted away from the origin, as well as an increase in \i for the first bands as 
x increases. This explains the increase of particle production as x shifts away from zero in 
Fig. 13. It is possible to compute the position of these bands based on (3.58) with (3.44), 
by determining the instability boundaries of (fto from the corresponding ones of Ak, but the 
approximation of small qi becomes unreliable as x and (fto grow. Thus, our analytic estimates 
provide a qualitative (and for small qi quantitative) explanation of the instability charts. 

The appearance of additional instability bands can be understood as well: by increasing 
x, we can no longer neglect the presence of q2, leading to a doubling of the period T and 
correspondingly the number of instability bands. 

The observed growth of /i with increasing x in the first bands leads to the dominance of 
the first one, and to a lesser degree the second one, during the late stages of preheating. To 
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Figure 15. The first instability-band of the k = 0-mode in a two dimensional held space, JV = 2, for 
different ESP distances from the origing, plotted over the initial values of the inflatons <f>o. 
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see this more clearly, it is instructive to associate a time t with each (fro according to (3.9), 
that is to identify the average oscillation amplitude with (fro so that 



As evident from the lower panel in Fig. 14, the zero mode spends the longest time in the 
first band. Thus, the first band can dominate preheating even though higher bands might 
be broader and have a larger characteristic exponent. From this plot we can also deduce the 
time at which preheating ends, since it corresponds to leaving the first resonance band once 
(fro becomes too small. 

We also observe a shift of instability bands to lower values of (fro as x is increased, 
resulting in an increased efficiency of preheating. We expect strong preheating should an 
instability band shift all the way down to (fro = 0. To see at which ESP positions such a 
shift takes place, we compute the first instability band in a two dimensional field space while 
increasing x from x = 0.0013 M p to x = 0.0036 M p in 16 increments, see Fig. 15. Evidently, 
the distance of the first band from (fr = is minimal for certain values of x*, leading to a 
prolonged (and thus strong) phase of particle production. As x is increased further, the band 
moves away from (fro = again, with ever decreasing \x up until it vanishes; finally, the next 
band repeats this sequence of events. 

Focussing on the first instability band, we see that it is closest to (fro = for x* « 
0.00141 M p , in accord with our analytic prediction in (3.60) for n = 2. 

However, we observed in Fig. 13 the absence of spikes for n > 4, opposite to the 
prediction of (3.60). To investigate the origin for their absence, we plot the instability bands 
corresponding to n = 4 and n = 5 in the lower two panels of Fig. 15: even though a given 
instability band approaches the origin as before, we observe an increased suppression of fi 
(and thus preheating) for small values of (fro that are particularly relevant during the late 
stages of the resonance phase. This explains the absence of spikes in Fig. 13. Similarly, no 
instability bands are present for n = 1, that is x\ = 7.071 x 10" 4 M p , and in case of three 
inflatons for n > 4. 

All of our analytic estimates were based on simple mass ratios, so that the period T is 
well defined. However, realistic models are expected to have more complicated mass spectra, 
such as the Marcenko-Pastur distribution in iV-nation [82]. 

How do our estimates fare with more complicated mass ratios? In Fig. 13, we observed 
the same qualitative behaviour, i.e. the presence of spikes, in line with our analytic expecta- 
tion if T is not identified with the period, but a characteristic time-scale of the oscillators. 
Thus, we still expect to see instability bands. To test this hypothesis, we computed the char- 
acteristic exponent in a two field model for two different mass ratios, mil 'mi = 1.5 (simple) 
and m2/m\ = y/2 (incommensurate) while keeping the overall scale unchanged, see Fig. 16. 

For small x, we observe a close correspondence between bands, especially for small (fro- 
Hence, the late stages of preheating are nearly indistinguishable in both cases. However, 
this correspondence becomes worse as x is increased and we loose the ability to predict the 
location of instability bands. Nevertheless, the qualitative features (clear instability bands 
are present that shift to lower (fro as x is increased) remain and we expect little change in the 
overall effectiveness of the resonance phase. 




(3.62) 
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Figure 16. Instability-bands of the k = mode for Af — 2 and x — 0, 0014 Mp (left) as well as 
x = 0, 00211 Alp (right). Solid line: simple mass ratio m-ijmx — 1.5; dotted line: incommensurable 
mass ratio vn^j vn\ = 



4 Summary and Discussion: Is Preheating Efficient? 

What is meant by efficiency? A pragmatic definition considers whether or not a seizable 
fraction of the inflatons' energy can be converted rapidly (compared to the Hubble time) 
into other degrees of freedom. In the presence of several inflatons, parametric resonance 
of a preheat matter-field \ that becomes light at the VEV of the inflatons is generally 
strongly suppressed due to de-phasing of inflatons [15, 38-40]; such de-phasing leads to a 
large additional contribution to the effective mass of %, preventing resonances. 

The main goal of this work was to investigate if this suppression prevails if one or 
more preheat matter fields become light at locations (ESPs) different from the VEV of the 
inflatons; the presence of several ESPs is generic for moduli spaces in string theory and often 
unavoidable in concrete models of inflation, e.g. in trapped inflation [43-46]. 

In Sec. 3 we observed and explained two main effects that can compensate for the 
suppression due to de-phasing, 

1. particle production during the the first in- fall (ESPs have to be densely distributed), 

2. particle production at certain ESPs (special locations) during a prolonged (narrow) 
resonance regime. 

We shall now summarize and discuss our main findings, before turning to lattice simulations 
to investigate the effects of back-reaction and double check our analytic results. 

4.1 The First In-fall 

During the first in-fall, particles are produced due to the violation of the adiabaticity condi- 
tion r\ < 1 with 7] in (2.14) at all ESPs that are within a certain distance of the trajectory 
in field space, roughly \(p — 0esp\ % fJ-c where fi c is the critical impact parameter in (3.24). 
We computed analytically the co-moving occupation number in the zero mode in (3.35) as 
well as the total number density in (3.40). We compared these results with numerics (no 
back- reaction) , see Fig. 11, and we found good agreement. Based on the analytic results, 
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we deduced that the suppression due to de-phasing can be more than compensated if ESPs 
are distributed densely and several inflatons are present, see Fig. 10 where we plot the total 
number density in the zero mode which we compare to the canonical single field preheating 
result (g = 5 x 10~ 4 and inflaton masses close to the COBE mass; for simplicity, we took the 
same coupling g between all pairs of \i an d </>i)- As a rule of thumb, the inter ESP distance 
should be below x < 0.005M p for J\f ~ 40: including more inflatons first increases and then 
decreases the total particle number (too many or too few inflatons suppress n t Q tal ); further, 
denser ESP distributions lead to more particle production. We checked that these results are 
insensitive to the concrete distribution of ESPs (only the average inter ESP distance mat- 
ters) or initial conditions of the inflatons (deviations are primarily caused by chance close 
encounters with individual ESPs, that boost particle production briefly). 

A word of caution: once the particle number becomes large, back-reaction needs to 
be incorporated: it is possible for the inflatons to slow down considerably, extending the 
inflationary phase; this effect is explained in detail in [46] and therefore not repeated here. 
It should be noted that lattice simulations are not feasible during the first in-fall, due to 
the plethora of preheat matter fields that need to be considered. Thus, back-reaction can at 
most be incorporated via the WKB approximation, as in [46]. 

Nevertheless, we can make some important, general observations: firstly, if this phase 
is efficient, a seizable fraction of the energy in the inflatons is distributed over many distinct 
preheat matter fields, but the occupation number in each Xi is l° w - Secondly, preheating can 
not be completed in this regime: once the speed along the trajectory falls below v ~ gx 2 , 
less than one ESP is (on average) close enough to allow for violation of the adiabaticity 
condition. Thus, at the very least energy of order v 2 ~ g 2 x 4 remains in the inflaton sector 
after the first in-flall. This energy needs to be transferred to other fields subsequently. A slow 
decay boosts the relative importance of these decay products, while a rapid decay, e.g. by 
resonance, retains the initial energy density ratio. 

4.2 The Resonance Phase 

After the first in-fall, a transitory phase with large (compared to the inter ESP distance x) 
amplitude oscillations follows, during which random bursts of particle production can occur 
whenever the trajectory comes close to an ESP, see Fig. 7. These intermittent bursts are usu- 
ally not sufficient to increase the particle number in preheat matter fields considerably. Once 
the oscillation amplitude falls below the inter ESP distance, a distinction into adiabatic/non- 
adiabatic regimes ceases to be meaningful, indicating that broad/stochastic resonances are 
impossible. Surprisingly, we found that a type of narrow resonance can be efficient in creating 
many particles in the zero mode if ESPs are at certain distances d n = x n ^/M /2 (with x n in 
(3.59)) from the VEV of the inflatons, see Sec. 3.4.1. These distances depend on the number 
of inflatons, their masses and the coupling(s) to the preheat matter fields. Further, a period 
T entered, which is only well defined for simple mass ratios (for incommensurate masses, it 
needs to be replaced by a characteristic time-scale of the oscillations). 

The presence of resonances can be understood intuitively by computing the character- 
istic exponent \i (/^ > indicates an exponential increase of no) for particle production at 
a single ESP in a static universe, see Fig. 14-16, where we plot u over the initial amplitude 
of inflatons 4>o (4>o is defined such that the total amplitude is A = (fo\0\f). As the location 
of the ESP moves away from the origin, instability bands move towards lower values of (po 
(whenever a band touches <po = 0, the ESP is at one of the special distances d n ). In an 
expanding universe, A decreases in time, scanning different bands in the process; naturally, 
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if a band reaches all the way down to </?o = 0, resonant particle production can operate for a 
long time and preheating is efficient. However, not all d n lead to efficient preheating, since 
the instability bands can become exceedingly flat (see the discussion at the end of Sec. 3.4.2); 
as a rule of thumb, only low n lead to strong resonances. 

The resonance phase ends once the oscillation amplitude becomes so small that the first 
resonance band is left 9 , see Fig. 14. Thus, to compute the time at which the plateau in Fig. 7 
is reached, one needs to identify first the ESP at which particle production is most efficient, 
compute the instability chart as in Fig. 14, read off the lower boundary of the first band, and 
finally plug this (fo into (3.62). 

A few comments are in order: particle production in this regime is dominated by a few 
ESPs that just happen to be at the right positions. If the average inter ESP distance in some 
random distribution is of order x\ or smaller, we expect that a few of them are at "right" 
distances. For the concrete case studies in Fig. 13, we find x\ ~ 1CP 3 , similar to the inter 
ESP distances at which the first in-fall is also efficient; as a consequence, we expect to see 
two distinct preheat matter classes: the first class consists of many different species with low 
occupation numbers, while the second one encompasses few preheat matter fields with high 
occupation numbers. 

In trapped inflation, ESP distributions need to be dense in order to cause any seizable 
back-reaction during inflation. As a consequence, the first in-fall is smoothly connected to the 
inflationary regime and we expect preheating in such models to be of the type examined in 
this article, including an efficient resonance phase. To decide if the inflatons decay completely 
after trapped inflation via resonances, one needs to perform lattice simulations; based on what 
we learned, we conclude that one does not need to consider all ESPs, but only those near 
the critical distances d n , which greatly reduces computational costs. We will not follow this 
line of research here: using a concrete implementation of trapped inflation in string theory is 
needed to include the necessary decay of preheat matter fields, which goes beyond the scope 
of the current paper. 

Nevertheless, we wish to put the above predictions with regard to the resonance phase 
to the test by means of lattice simulations in a few simple scenarios in Sec. 5, primarily to 
make sure that resonances persevere in the presence of back-reaction, but also to check some 
of the approximations that went into the analytic results. 

5 Lattice Simulations 

To solve the coupled system of the full field equations in an expanding background, we use 
the same modified LATTICEEASY code as in [40] (see Appendix A of [40] for details on the 
code implementation as well as initial conditions for inhomogeneous modes). As mentioned 
in Sec. 2.2.2, we use a coupling of gi = g = 5 x 10 -4 and start the matter field(s) in their 
vacuum state. Since we will not be studying any tachyonic preheating scenarios as were 
studied in [40], we neglect the self interaction for the matter fields, Xxf- With this tool at 
hand, we wish to discuss two setups: 

1. Particle production at a single ESP, with varying number of inflatons, initial conditions 
and masses, to test the predictions of Sec. 3.4, particularly the viability of a prolonged 
narrow resonance regime as an efficient preheating mechanism. 

9 If the first band reaches all the way down to ipo, we still have n(tpo) — > smoothly, and particle production 
still becomes inefficient at some point. 
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Figure 17. The ratio of the energy in the \ field compared to the energy in <f> for Af = 2. The color 
coding shows x = 0.00106 M p (red), x = 0.00141 M p (orange), x = 0.00177 M p (blue), x = 0.00212 M p 
(brown), and x = 0.00247 M p (black). 
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Figure 18. The ratio of the energy in the \ field compared to the energy in 4> for J\f — 2 and a 
single ESP at the origin, but a bare mass of the matter field, = g 2 x 2 /2. The color coding shows 

x = 0.00106 M p (red), x = 0.00141 M p (orange), x = 0.00177 M p (blue), x = 0.00212 M p (brown), 
and x = 0.00247 M p (black). 

2. Particle production at several ESPs, Gaussianly distributed near the origin (see Sec. 2.2.2), 
to test if preheating is indeed dominated by a few matter fields. 

5.1 Particle Production at a Single ESP 

The first question we would like to address is whether or not the prolonged period of narrow 
resonance during small amplitude inflaton oscillations is sufficient to transfer 0(1) of the 
inflatons energy into a matter field. We choose the same setup as in Sec. 3.4, varying the 
ESP position by changing x over an interval that covers the spikes in Fig. 13. To test the 
efficiency, we compute the ratio of the energy in the matter field p x over the total energy in 
the inflaton sector p v . 

All simulations in this section use a 128 3 lattice, each side has initial length Lq = 5m _1 , 
and we run the simulations until t = 400m~ 1 unless otherwise noted. 

In Fig. 17, we plot the energy ratio for two inflatons with m^/mx = 3/2 , mi = m = 
10~ 6 Mp, we use the initial field values denned by Eq. (2.8), zero initial speed and a single 
ESP at ES p = (x,x)/2. We vary x to cover the region around the first two peaks at 
x w 0.00141 M P and x 0.00212 M P in Fig. 13 top left panel (see equation (3.59) for the 
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Figure 19. The ratio of the energy in the x field compared to the energy in <j) for JV = 3. The color 
coding shows x = 0.00115 M p (red), x = 0.00130 M p (orange), x = 0.00144M p (blue), x = 0.00159 M p 
(brown), and x = 0.00173 M p (black). 

analytic value). We observe that as long as x is in the vicinity of the peaks, resonances are 
indeed sufficient to raise the ratio p x / 1 p v to order one, at which point the fields fragment and 
subsequently enter a turbulent regime. The exact value of x has little effect onto preheating, 
as long as it is in the region of the peaks. This is expected, since altering x changes primarily 
the maximal duration of the narrow resonance regime, without significantly changing the 
generalized Floquet index. Hence, as long as resonances last long enough to bring this ratio 
to order one, preheating operates similarly. However, once the ESP is too far away, i.e. x is 
to the right of the second peak in Fig. 13 top left panel, resonances shut of before preheating 
is completed (black line in Fig. 17). It should be noted that an ESP at the origin can still 
lead to effective preheating for J\f = 2, see i.e. Fig. 1A in [40] or Fig. 4, such that prolonged 
narrow resonances may not be crucial if only two inflatons are present. 

It is instructive to compare the observed resonances to the ones with a single ESP at the 
origin, but allowing for a non-zero mass of the matter field w^g- = g 2 x 2 /2, see Fig. 18; the 
initial stages of parametric resonance are also described by the analytic result in (3.59), since 
we took the limit <3?o — > and consequently ignored q\ and ■ As expected, we observe 
that particle production becomes less efficient as the mass of the x- field grows, and ceases to 
be important if 



In our case resonances become more and more suppressed once the corresponding x is to 
the right of the first peak in Fig. 13 (top left panel); for x > 0.0025 M p (x to the right of 
the second peak) resonances are hardly present at all. Interestingly, preheating is efficient 
in Fig. 17 with the corresponding choice of ESP position: thus, to explain p x /p v ~ 0(1) 
in Fig. 17, the inclusion of q£\ q£ and/or backreaction are important, while our simple 
analytic treatment remains only qualitatively valid. We conclude that lattice simulations are 
crucial: the presence of resonances in perturbative analytic calculations are necessary, but 
not sufficient, indicators to ascertain the efficiency of preheating. 

Let us turn our attention to the J\f = 3 case, Fig. 19, for which parametric resonance 
with an ESP at the origin is strongly suppressed due to dephasing, see Fig. 1A in [40] or 
Fig. 4. Choosing mj/mi = 5/4, m^/mi = 3/2, mi = m = 10~ 6 Mp and (pi n i = (l,l,l)x/2, 



m xcff > f ew x 9 2x i/2 = few x 




(5.1) 
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Figure 20. The ratio of the energy in the \ ne ld compared to the energy in <f> for Af = 2. 
The color coding shows x = (0.00141, 0)/y/2M p (red), x = (0, 0.00141)/y^M p (blue), x = 
(0.00141, 0.00141)/2M P (orange). 



1 F 




50 100 150 200 250 300 
Rescaled Time (m~ l ) 



Figure 21. The ratio of the energy in the \ held compared to the energy in <f> for N = 2. The color 
coding shows case (a) where m^lmy = 3/2 for x = 0.00141 M p (red) and x = 0.001 M p (orange) and 
case (b) where m 2 /™i = V% for x — 0.00141 M p (green) and x = 0.001 M p (blue). 

corresponding the top right panel in Fig. 13, we observe again efficient, practically indis- 
tinguishable, preheating for x near the peaks in Fig. 13. Once the last peak is traversed 
(x > 0.00173 M p , the black line in Fig. 19), narrow resonances weaken, shut off earlier and 
preheating becomes inefficient. 

So far we kept the ESP at </? ESP = (x, . . . , x)/2, just changing the value of x. Based on 
our analytic derivation leading to the prediction of the peaks' positions in (3.59), we do not 
expect a strong dependence on the angular position of the ESP if the oscillation amplitude 
and Xi are small enough to ignore (defined in (3.44)). We can confirm this expectation 
in Fig. 20, where we keep the distance of the ESP to the origin fixed (at the first peak), but 
change its angular position. While some angular dependence is evident regarding the onset of 
resonance, the final amplitude of p x /Pip remains the same. If the ESP is pushed out further, 
we expect to see a stronger breaking of spherical symmetry, see Sec. 5.2. 

The analytic prediction of the peaks position in (3.59) was based on the assumption of 
commensurate masses. However, as we argued in Sec. 3.4.1 and saw i.e. in the lower panels 
of Fig. 13, the qualitative features of the resonance phase remain valid for incommensurate 
masses; this feature remains true on the lattice, see Fig. 21, where efficient preheating via 
narrow resonance is still present. 
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Figure 22. The location of ESPS in ten different simulations containing ten ESPs each. The left 
hand panel shows the most-significant ESPs that contribute to 60% of the total energy of the x-fields 
(blue squares) in each simulation compared to the less energetic ESPs (red circles). The right hand 
panel shows the most-significant ESPs that contribute to 90% of the total energy of the x-fields in 
each simulation (blue squares) compared to the less energetic ESPs (red circles). For comparison, 
the two circles represent the locus of points that correspond to x = 0.00141 (inner) and x — 0.00212 
(outer). 

We now want to look at resonance when more than one ESP is present in the simulation. 
First, we examine the case of N = 2 inflatons. Due to time limitations on our simulations, 
we consider independent runs with 10-ESPs per simulation. The locations of the ESP's are 
drawn from a Gaussian distribution with variance a = 1.875 x 10 _4 M^, so that there 
are a number of ESP's near the peaks of Fig. 13. Again we chose the ratio of masses as 
■m-zjm\ = 3/2, mi = m = 10~ 6 Mp, with initial field values defined by Eq. (2.8) and zero 
initial speed. 

In Fig. 22 the position of ESPs in ten of these simulations is plotted accumulative on the 
same panels. ESPs at which matter fields are predominantly produced are indicated by blue 
squares (making up 60% and 90% of the total energy in the left and right panels respectively). 
The two circles indicate the position of the first and second peak in Fig. 13 (top left). Firstly, 
we observe that ESPs further out dominate, particularly near the second circle, while ESPs 
near the origin are unimportant. This is in line with our analytic and numeric conclusions 
in previous sections. Secondly, we observe the breaking of spherical symmetry, particularly 
in the right panel: ESPs can be closer in the 932 direction (the heavier field) than in the tpx 
direction, while remaining important for preheating. This breaking of spherical symmetry is 
caused by q!p in (3.59), which we ignored in the analytic discussion. As we saw by comparing 

(i) (i) 

Fig. 17 to Fig. 18, the presence of q\ and q\ is important to attain a large energy fraction. 
It is thus not surprising that ESPs dominating preheating are the ones that show a breaking 
of spherical symmetry. 

A similar breaking of spherical symmetry should be present for M > 3. However, even 
the N = 3 case is quite time consuming to treat on a lattice, since it requires a great number 
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Figure 23. A histogram showing the number of ESPS compared to the distance from the origin in 
field space for ten simulations containing ten ESPs each. Each ESP is located at </?esp = (x,x,x)/2, 
where x is drawn from a Gaussian distribution. Three most-significant ESPs that contribute to the 
energy of any given simulation (blue) are compared to the total number of ESPS at that distance 
(red). 



of ESPs to yield an adequate filling of field space. Hence, we only investigate a Gaussian 
distribution of ESPs along a prescribed axis (again accumulative, with 10 ESPs per run): in 
Fig. 23 we chose </?esp = (x,x,x)/2, with x Gaussianly distributed. The histogram shows 
the ESPs that contribute the most to preheating (60%, violet), while the less important ones 
are coloured pink. We observe again that preheating is not dominated by ESPs at the origin, 
but at distances bigger than the fist peak in Fig. 13 (top right), in line with the discussions 
in previous sections. It is tempting to identify peaks in Fig. 23 with the ones in Fig. 13, but 
the total number of ESPs in our simulations, and thus the statistical significance, is too low 
for such a detailed comparison. 

6 Conclusion 

Motivated by the presence of extra species points (ESP) in many models of inflation in string 
theory (i.e. trapped inflation), we investigated preheating of several preheat matter fields 
that become light near, but not at, the vacuum expectation value (VEV) of the inflatons. 
In the presence of several inflaton fields, dephasing tends to suppress resonances, which can 
however be compensated by two new, hitherto overlooked, effects: first, if ESP distributions 
are dense, particles can be produced at many ESPs via broad resonance once slow roll ends, 
since the inflaton speed is not slow roll suppressed any more and the adiabaticity condition is 
violated for many preheat fields. We compute analytically the resulting occupation number, 
paying particular attention to the dependence of preheating on the number of inflatons and 
the density of ESPs. Preheating during this first in-fall, albeit potentially sufficient to convert 
the majority of the inflatons energy to preheat matter fields, is however never complete. 

After an intermediate phase of large amplitude oscillations during which additional 
particle production may occur, a prolonged phase of narrow resonance generically leads to the 
complete decay of inflatons, if the location of an ESP is at one of several "right" distances from 
the VEV. We compute these locations analytically, by means of a simple stability analysis in 
a static universe, providing a thorough, yet intuitive understanding of the resonance regime. 
We test our analytic predictions numerically by direct integration of the equations of motion 
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(no backreaction) and complementary lattice simulations and find good agreement with the 
analytics. 

Our study shows that the understanding of preheating after multi-field inflation is far 
from complete: we found new, efficient mechanisms of preheating that can be missed easily 
if the model is over-simplified. For example, while the presence of ESPs in field space is 
often incorporated during inflation, they are usually ignored or lumped together into a single, 
effective preheat matter field that becomes light at the VEV of the inflatons, which completely 
eliminates the preheating mechanism discussed in this paper. Further, combining inflatons 
into an effective single field, albeit extremely useful and predictive during inflation, alters 
the nature of preheating qualitatively and is therefore never a prudent approach towards the 
end of inflation. 

One of the main, remaining challenges for inflationary model building in string theory 
is the identification of all relevant degrees of freedom that the inflaton(s) couple to once 
inflation terminates. Our work is a step in this direction by incorporating the presence 
of ESPs, describing fields that are present in almost all moduli spaces of string theory. 
Consequently, the next step towards a theory of reheating after inflation in string theory 
is the incorporation of subsequent decay products and the computation of thermalization, 
which we leave for future studies. 
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